Diffusive and ballistic motion in superconducting hybrid structures. 
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We examine transport properties of superconducting hybrid mesoscopic structures, in both the 
diffusive and baUistic regimes. For diffusive structures, analytic results from quasi-classical theory 
are compared with predictions from numerical, multiple-scattering calculations performed on small 
structures. For all structures, the two methods yield comparable results and in some cases, quan- 
titative agreement is obtained. These results not only demonstrate that quasi-classical theory can 
yield the ensemble averaged conductance < G > of small structures of dimensions of order 10- 20 
Fermi wavelengths, but also establish that numerical scattering calculations on such small structures 
can yield results for < G > which are characteristic of much larger systems. Having compared the 
two approaches, we extend the multiple-scattering analysis to the ballistic limit, where the sample 
dimensions become smaller than the elastic mean free path and demonstrate that the properties of 
certain Andreev interferometers are unchanged in the clean limit. 



Pacs numbers: 72.10B, 73.40G, 74.50 
I. INTRODUCTION 

When a normal metal (N) makes contact with a super- 
conductor (S), classical tunnelling theory predicts that as 
a consequence of the existence of a superconducting en- 
ergy gap, the low-temperature sub-gap conductance will 
vanish. Andreev scattering |^ provides an alternative 
mechanism for charge transport through such an N-S 
junction and for an ideal interface Q], leads to a zero- 
voltage conductance which is almost twice that of the 
normal state. In the presence of an insulating layer (I), 
Andreev scattering becomes less effective and in the ab- 
sence of disorder, the sub-gap conductance of an N-I-S 
junction is predicted to be depressed compared with that 
of the normal state. In contrast, at low enough temper- 
atures, experiments on N-I-S junctions Q reveal the ex- 
istence of a zero- voltage conductance peak, with a value 
comparable to the conductance in the normal state. This 
effect is due to the interplay between disorder-induced 
scattering and tunnelling and has been observed in ex- 
periments involving quantum wells with superconduct- 
ing electrodes Q, superconductor- normal metal micro- 
junctions H, and superconductor-2DEG-superconductor 
structures^. These have been interpreted using a num- 
ber of theoretical approaches, including quasi-classical 
Green function techniques |@]- ||l^, tunnelling Hamil- 
tonian methods ||l^, [Q, multiple-scattering techniques 
||l3, and random matrix theory |jl|], 

Recently, following a number of theoretical papers on 
disordered transport in the presence of two superconduct- 
ing contacts |^^- |2^, several new experiments aimed at 
probing the phase-coherent nature of Andreev scattering 
have been carried out. These involve a normal metal in 
contact with superconductors 5*1 and ^2, with order pa- 
rameter phases <j)i, 4>2, whose order parameter phase dif- 
ference (j) = — (j)2 can be controlled by external means. 
For a diffusive system of size L, with diffusion coefficient 
D, early theoretical work predicted that in the 



high temperature limit T > T*, where ksT* = hD/L'^, 
the ensemble averaged conductance < G > should be 
a periodic function of (j), with fundamental period tt, 
whereas in the low temperature regime T < T*, it was 
predicted |22), |^ that the fundamental periodicity of 
< G > would be 2tt. This prediction of a 27r periodicity 
is a common feature of all recent theories of < G > , 
|l|, H, H and of recent experiments 1^- 

Despite much progress, many details of such Andreev 
interferometers remain to be understood. In particular, 
there exists no general theory of the amplitude of oscilla- 
tion, the nature of the zero-phase extremum and the har- 
monic content of the conductance-phase characteristic. 
Experiments on various geometries, in different transport 
regimes have yielded amplitudes of oscillation which dif- 
fer by many orders of magnitude. It has also been ques- 
tioned whether analytical theories that mainly apply in 
the diffusive limit can be used in the quasi-ballistic case 

One difficulty in establishing a general theory is that 
most theoretical papers are based on a single technique, 
with little detailed comparison with the results of other 
approaches. For example, distinct analytical theories 
exist for the ballistic, diffusive and strongly disordered 
regimes, but there is no analytical theory capable of de- 
scribing the cross-over between them. In contrast, exact 
numerical solutions of the Bogoliubov - de Gennes equa- 
tion [^2| can easily cross from one regime to another, 
but are limited to system sizes of order a hundred Fermi 
wavelengths. 

In this paper, we undertake a detailed comparison be- 
tween an exact numerical, multiple-scattering technique 
and quasi-classical theory. The former was first used to 
solve the Bogoliubov - de Gennes equation for disordered, 
one-dimensional systems [ p3[ and soon after generalized 
to higher dimensions |^^, while the latter, developed 
in the context of nonequilibrium superconductivity [|35[|- 
p8|, when supplemented by the appropriate boundary 
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conditions [g9| , |Q , has recently yielded a variety of re- 
sults for ensemble averaged conductances in the diffu- 
sive limit 0- Once agreement between the two ap- 
proaches is obtained in the diffusive limit, we then use the 
numerics to follow the crossover to the ballistic regime, 
where the electron mean free path becomes comparable 
with the system size. 

The comparison will be carried out by examining two 
canonical examples of mesoscopic superconducting hy- 
brid structures, namely an N-I-S junction and an An- 
dreev interferometer. In section II we briefly recall some 
results of the quasi-classical Green function approach for 
these structures and in section III highlight the main fea- 
tures of the numerical scattering approach. Whereas a 
quasi-classical approach yields a one-dimensional theory 
whose results depend only on the topology of the struc- 
ture, the numerical calculations, in common with real 
experiments, require the specification of a suitable geom- 
etry. In section IV we present a detailed discussion of the 
geometry and numerical parameters needed to reproduce 
the results of quasi-classical theory. Having established 
agreement between the two approaches we then depart 
from the diffusive limit and investigate the cross-over to 
ballistic transport, which is particularly relevant to the 
experiments of reference [pO[. 



II. RESULTS FROM QUASI-CLASSICAL 
THEORY. 



In this section, we highlight some predictions of the 
quasi-classical approach of references [Q- ||T^. These the- 
ories focus on the ensemble averaged conductance < G > 
and ignore the weak localisation contribution discussed 
in |2^. The latter does not scale with the system 
size and in systems with a conductance much larger than 
2e^/h, can be neglected. The results of references 0- 
1^ are obtained by solving the following equation for 
the quasi-classical Green function g in the diffusive limit 
Ddr{gdrg) + iE\fz, g] = 0, where E is the quasi-particle 
excitation energy |41[| . In what follows, we shall consider 
only the E — limit of this equation, which applies to a 
diffusive system, whose length L is assumed to satisfy the 
inequalities L <^ ^, where = D/Aq and — D/E, 
in the superconducting and normal region respectively. 
Here D is the diffusion coefficient, Aq the energy gap 
and E — eV, where V is the applied voltage. Further- 
more all temperatures and voltages are assumed to be 
much smaller than Aq. For convenience in what follows, 
we shall also restrict the analysis to zero temperature. 

In units of 2e^ /ft, theory predicts jT|], that the total 
conductance G of the structure shown in figure |f|(a) may 
be computed from the equation 



I 
G 



I 



1 



Gd. 



iff 



Gtun sin 9 



(1) 



where is a solution of the transcendental equation 

Gdiff 9 = Gtun COS 9 



(2) 



and Gtun j Gdi / / are the respective conductances of the 
tunnel junction and diffusive region in isolation. 

There are two obvious limits to take. The first is where 
Gtun /Gdiff for which 9 « Gtun /Gdiff « sin 6* and 
hence 
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The second limit is Gtun/Gdif f 
7r/2. This yields 
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Eq.(^ has been also directly obtained in the tunnel- 
ing Hamiltonian limit p^ . Furthermore, as emphasized 
in Ref. the change of the power in the dependence 
on the tunnel junction conductance Gtun, as described 
by the equations (H) and , reflects the combined effect 
of the Andreev scattering at the N-S interface and the 
interference effects in the mesoscopic phase coherent dis- 
ordered normal region. In the regime described by eq. (Q) 
and when Gtun *S 1, the conductance of the N-S junc- 
tion appears to be the same as in the normal state and 
the total resistance is obtained by simply adding up in 
series the resistances of the diffusive region and tunnel 
junction. Equations (|l|) and (H) may be solved numer- 
ically to yield the resistance Rtot — 1/G of the system 
as a function of Gdiff and Gtun- The result of such an 
exercise for a fixed value Gdiff = 1-6 and a variable Gtun 
is shown by the curve (a) of figure E. Also plotted are the 
limits given by equations (^) and (|) represented by the 
dashed fines (e) and (d) . The curve (b) shows the results 
of the numerical simulation described in section IV. 

As a second example, quasi-classical theory predicts 
| pT[ that the conductance G of the structure depicted in 
figure ||(a) may be computed from the equation 
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where Gi is the conductance of the tunnel junction (1), 
G2 is the conductance of the tunnel junction (2) and (j) 
is the phase difference between the two superconductors. 
In the limit Gi ^ G2, this simplifies to the expression 
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whereas if G2 ^ Gi then 
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Figure ^ shows four plots of equation (||) for each com- 
bination of Gi = 0.2, 2.0 and G2 = 0.2, 2.0. One can see 
that for Gi — G2 there is a zero-phase minimum, as there 
is for Gi ^ G2. When Gi ^ G2 however, one finds a 
zero-phase maximum. In all cases the conductance van- 
ishes when the phase difference between the supercon- 
ductors is TT. Corresponding analytical results have been 
obtained for ballistic interferometers in one |Q and 
two [H dimensions. 



III. A SCATTERING APPROACH TO 
TRANSPORT IN MESOSCOPIC 
SUPERCONDUCTORS. 

During the past six years ||3^- jsj] we have developed 
numerical codes capable of yielding exact solutions of the 
Bogoliubov - de Gennes equation for disordered struc- 
tures in arbitrary dimensions. Currently there are two 
independent sets of codes available at Lancaster; one of 
these is based on a transfer matrix approach and the 
other is based on a recursive Green's function method. 
Typically these are used as independent cross-checks and 
therefore there can be no doubt about the accuracy of 
the results obtained for a given structure. For a two- 
dimensional system of width less than a few hundred 
Fermi wavelengths, or for a three-dimensional system of 
width less than a few tens of Fermi wavelengths, the 
problem of computing dc transport properties of a phase- 
coherent system described by mean-field BCS theory is 
therefore no longer an issue. Just as the appearance 
of pocket calculators rendered approximate methods for 
computing elementary functions redundant, the existence 
of these codes has, for several years, allowed transport 
properties of small structures to be calculated without 
further approximation. For larger systems, the key is- 
sue is how to extrapolate the results of such calculations 
to larger numbers of channels. By making contact with 
quasi-classical theory, the results which follow establish 
that for many systems, the ensemble averaged conduc- 
tance obtained from small systems exhibits the essential 
features of much larger structures. 

The numerical codes evaluate multi-channel scattering 
formulae for the dc electrical conductance [|l5| js^ and 
more recently thermoelectric properties ||45| . In what fol- 
lows, we focus on hybrid structures connected to normal 
external reservoirs only |46|. For a structure connected 
to two such reservoirs, the zero-temperature, zero-bias 
electrical conductance can be written [|l5|, (in units 
of 2e^/h), 



G = To+Ta 



2{RaR'^ - TgT^) 

Ra + R'+Ta + T' 



(8) 



In this expression, i?o,7o {Ra,Ta) are the coefficients 
for normal (Andreev) reflection and transmission for 



zero-energy quasi-particles from reservoir 1 , while R'q , Tq 
(i?^,T^) are the corresponding coefficients for quasi- 
particles from reservoir 2. If each of the external leads 
connecting the reservoirs to the scatterer contains N 
open channels, these satisfy Ro+To + Ra+Ta — Rq+Tq + 
R'^ + T'^ = N and Tq + = T^ + T^- Furthermore, in the 
absence of a magnetic field, all reflection coefficients are 
even functions of <j), while the transmission coefficients 
satisfy T^((/.) = To{^c^), T'M = Ta{-cp). Consequently 
on quite general grounds, in the absence of a magnetic 
field, G is predicted to be an even function of 0. 

In the absence of quasi-particle transmission between 
the two external probes, equation (H) reduces to 
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where 2Ra {2R'^) are left (right) boundary conductances, 
introduced by Blonder, Tinkham and Klapwijk j^]. The 
Lancaster codes yield all possible scattering coefficients, 
but in what follows, to compare results with the theory 
leading to eqs.([^) and (j|), we shall analyze structures 
with no transmission and restrict attention to Ra only. 

The numerical codes compute the scattering coeffi- 
cients of a tight-binding lattice, described by a Bogoli- 
ubov - de Gennes (BdG) Hamiltonian of the form 
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A* 
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(10) 



If an index n is used to label sites on the lattice and 
any internal spin-degrees of freedom, then Ho is of the 
form 



Ho = ^e„|n)(n| + ^ Vn,m\n){m\. 



(11) 



(n.m) 



To model a given physical structure, it is necessary to 
specify certain phenomenological parameters which cap- 
ture the essential physics. As an example, in the absence 
of spin-orbit scattering, spin degrees of freedom can be ig- 
nored and in the absence of a magnetic field, one chooses 
Ki,m — —7 for nearest neighbour pairs {n,m). If (n,m) 
are not nearest neighbours, then Vn^m = 0. In a region 
free from disorder, the diagonal elements e„ are set equal 
to a constant eo: whereas in a disordered region, e„ is a 
random number uniformly distributed between eo — W 
and Co + W. In the presence of spin singlet, local s-wave 
pairing, A is a diagonal order parameter matrix with el- 
ements A„. In a normal region, A„ = 0, whereas in a 
clean superconducting region, |A„| is set to a constant 
value Aq. The phase of A„ is chosen to equal a value 
assigned to the superconducting region to which site n 
belongs. In what follows, the energy scale will be fixed 
by making the choice 7=1. 

Of course, the above parameters are not directly ac- 
cessible experimentally and are not an explicit feature 
of quasi-classical theory. Therefore when making com- 
parisons, some effort is needed to map one analysis onto 
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another. In d-dimensions, for a clean system on a square 
or cubic lattice, the chemical potential relative to the 
band bottom is ^ = eq + 2(1^, the band width is and 
the effective mass for excitations near the band bottom 
is m* = ?i^/(27a^), where a is the lattice constant. A 
key parameter in the problem is the dimensionless ratio 
A = Ao//i, which takes a value 10~^ for conventional 
low Tc superconductors such as Niobium, but can be 
as large as 0.1 for high-temperature superconductors, or 
for a 2DEG in contact with a conventional superconduc- 
tor. Andreev's approximation, which underpins many 
analytic theories, including quasi-classical and random 
matrix descriptions, is valid only when this parameter is 
much less than unity. 

Other parameters which are needed when making com- 
parisons are the Thouless energy E*, which for a dif- 
fusive structure of width M, length Ldiff and normal- 
state conductance G, is given by E* = hD/L'^^JJ^ — 
{h/2e^)G/{n{0)LdiffM'^-^), where n(0) is the density of 
states per site. A second parameter is the normal-state, 
elastic mean- free path I, which for a diffusive sample con- 
nected to external lead with N open channels, is given 
by G = {2e^ /h)Nl/ Ldiff ■ Within a numerical simulation 
on a given geometry, once the model parameters W ^ 7, 
eo and Aq are chosen, the parameters G, Z, n(0) and E* 
are computed explicitly. 

IV. NUMERICAL RESULTS FOR A N-I-S 
STRUCTURE. 

In this section, we present a comparison between the 
predictions of quasi-classical theory and the above nu- 
merical scattering approach, for an N-I-S structure. Our 
aim is to highlight the steps required to obtain a suitable 
choice of parameters, from which a meaningful compar- 
ison can be made. In the literature, numerical results 
in two-dimensions have been obtained by first solving 
for the scattering matrix of a normal diffusive structure, 
with or without a tunnel junction and then employing 
Andreev's approximation to model the Andreev scatter- 
ing induced by a nearby superconductor As noted 
above, this approximation requires that Aq be small com- 
pared with the Fermi energy and that there be no dis- 
order in the superconductor. Furthermore for a clean 
system, Andreev's approximation can yield incorrect re- 
sults, because even at a clean N-S interface, the approx- 
imation breaks down when scattering channels are 
almost closed. For these reasons a comparison with an 
exact solution of the Bogoliubov - de Gennes equation 
allows one to examine changes occurring away from the 
Andreev limit, in a region of parameter space which is 
more relevant to high-temperature superconductors. 

The system to be examined is shown in figure 0(b) 
and consists of a disordered region in contact with a 



tunnel junction, which is in turn adjacent to a super- 
conducting probe. The simulated structure consists of a 
two-dimensional tight-binding lattice of width M sites. 
The disordered region is of length Ldiff sites, the tunnel 
junction is Lt„„ sites long and the superconductor has a 
length Lsup- In units of 2e^//i, the conductance of a par- 
ticular realisation of the structure will be denoted G and 
the ensemble-averaged conductance will be written (G). 
In order to make a comparison with quasi-classical the- 
ory, it is necessary that the properties of the two resistive 
components and the superconducting probe be compati- 
ble with the assumptions made by the theory. 

To identify a suitable choice of parameters, consider 
first a normal diffusive portion of length Ldi / / and width 
M, connected to crystaUine leads. In units of 2e^//i, the 
conductance of a particular realisation of this structure 
will be denoted Gdiff and the ensemble-averaged con- 
ductance will be written {Gdiff)- The conductance of a 
diffusive material is inversely proportional to its length 
and therefore a plot of {Gdi ff)Ldi // as a function of Ldi / / 
will exhibit a plateau in the diffusive regime, with a mean 
free path given by / = {Gdif f)Ldif f /N . 

For a sample of width M = 10, figure || shows a plot 
of {Gdif f)Ldif f /N versus Ldiff - This structure has peri- 
odic boundary conditions in the direction tranverse to 
the current flow and the choice eo = 0.2 was made, 
which yields = 9. Results are shown for a disorder 
of = 1. The length Ldiff of the disordered region 
was incremented in steps of two sites from Ldi // = 2 to 
40. For each value of Ldiff, 2000 realizations of disor- 
der were chosen and the conductance Gdiff computed 
for each. Then the ensemble average {Gdiff) and the 
standard deviation SGdiff were calculated. 

Figure |^ shows that in the interval 20 < Ldiff < 40 
the system exhibits diffusive behaviour, with a mean free 
path of I « 4.6. For smaller values of Ldiff, the system is 
in the ballistic regime and for larger values, the onset of 
localisation causes the curve to fall. A diffusive system 
is one for which I <C Ldiff and I <^ M - Furthermore 
if weak localization corrections are to be neglected, we 
require Nl 3> Ldiff - For these reasons a judicious choice 
of length yielding a diffusive system whilst minimizing 
CPU time is Ldiff = 25. Such a system has an average 
conductance {Gdiff) = 1-6 

To compare with quasi-classical theory, a knowledge of 
the conductance of the isolated tunnel junction Gtun as a 
function of the barrier height £{, is also required. In what 
follows, we consider a clean tunnel junction of dimensions 
M — 10 and Ltun = 1, obtained by setting the diagonal 
elements e„ of all barrier-sites n equal to eo -I- Cfc. For an 
isolated barrier connected to crystalline leads of width 
10, figure ^ shows a plot of Gtun as a function of e^. To 
obtain this plot, the conductance Gtun is computed for 
1000 successive values of the barrier height eh in the range 
0.0 < e;, < 10.0. This choice of barrier heights yields a 
spread of barrier conductances in the convenient range 
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0.0 < Gtun < 9.0. 

Finally, before a comparison with theory can be made, 
the properties of the superconductor must be chosen 
such that there be no transmission through the super- 
conducting region and that Andreev's approximation of 
neglecting normal reflection at the N-S interface is valid. 
To avoid quasi-particle transmission, it is necessary to 
choose Lsup > where ^ = /^/Aq and to minimise nor- 
mal reflection it is necessary that ^ >> 1. For Aq = 0.05 
and eo = 0.2 the superconducting coherence length is 
^ = 76 and it is found that transmission is negligible for 
Lsup > 100. For the above choice of parameters, one finds 
for the normal and Andreev reflection and transmission 
coefficients: Ro = 0.06485, Ra = 8.84762, To = 0.08748, 
and Ta = 0.00005. 

It should be noted that the condition ^ >> 1 is not 
sufficient to completely exclude normal reflection at an 
N-S interface. It is also necessary that eo be chosen 
such that the number of open channels in the external 
leads is not sensitive to small changes in eo. This fea- 
ture is illustrated in figure 0, which shows as a function 
of eo, the conductance of a clean superconducting re- 
gion of width M = 50 and length Lsup — 5, attached 
to crystalline normal leads. Results are shown for 5 val- 
ues of Ao: Ao = 0, 10-2, 10-\ 0.3 and 0.5. For Ao = 0, 
the conductance is equal to the number of open chan- 
nels and changes by unity whenever an external quasi- 
particle channel closes. At these values of eo, switching 
on an infinitesimal Ao causes the conductance to decrease 
by unity. As shown in the figure, finite values of Ao 
smear the conductance steps and suppress the conduc- 
tance. Both of these features lie outside Andreev's ap- 
proximation. To achieve compatibility with the assump- 
tions of circuit theory, the choice eo = 0.2 was made, 
which places a system of width 10 between two conduc- 
tance steps and avoids the above sensitivity to changes 
in Ao. 

Having identified a choice of parameters which is com- 
patible with quasi-classical theory, numerical results for 
the combined structure of figure ^b can now be ob- 
tained. To summarize, this structure has the follow- 
ing properties: width M = 10, number of open chan- 
nels N = 9, band filling factor eo = 0.2 leading to a 
chemival potential ^ = 3.8, length of tunnel junction 
Ltun = 1, barrier heights 0.0 < ef, < 10.0, barrier con- 
ductances 0.0 < Gtun < 9.0, length of diffusive region 
Ldiff = 25, diffusive disorder width W = 1, conduc- 
tance of diffusive region (Gdiff) — 1-6, length of super- 
conductor Lsup — 100, superconducting coherence length 
^ = 76, superconducting order parameter Ao — 0.05, 
elastic mean free path I « 4.5. 

First consider the case of no barrier, where e^ — 0. 
In this case, quasi-classical theory insists that the con- 
ductance of a diffusive region in contact with a super- 
conductor should be identical with the normal-state con- 
ductance of the diffusive region. Figure p| shows plots 



of the mean conductance (G) = N — i?o(Ao) + i?a (Ao) 
as a function of Ao, for disordered regions of four dif- 
ferent lengths. In the normal state (Ao = 0) (G) 
reduces to T'o(O) = N — i?o(0) and in the presence 
of a sufficiently-long superconductor, to the BTK con- 
ductance 2i?Q(Ao). The left insert shows the quan- 
tity < 5 (iV - i?o(Ao) + i?a(Ao))/(To(0)) (ie the 
conductance divide by the normal state conductance). 
The right insert shows the root mean square deviation 
a = ((G(Ao) - {G{Aa)))y/^. These show that in the 
ballistic limit Ldiff = 0, the conductance rises to a value 
almost double that of the normal state, before decreasing 
with increasing Ao. In contrast, the mean conductance 
of a diffusive normal region is relatively insensitive to 
the onset of superconductivity, with the largest relative 
change corresponding to the largest value of Ldiff, ( ie 
the smallest value of the normal state conductance) . It 
should be noted however that even though the mean is 
insensitive to Ao, the fact that the rms deviation a is non- 
zero reveals that for individual samples, large changes of 
arbitrary sign can occur. 

Having examined a diffusive conductor with no bar- 
rier, we now turn to the case of finite ef,. Curve (b) of 
figure 1^ shows numerical results in the presence of a tun- 
nel barrier. For 50 equally-spaced barrier heights in the 
range 0.0 < e;, < 10.0, 500 realizations of disorder in the 
diffusive region were selected and the total conductance 
G computed for each realisation. The ensemble-averaged 
conductance (G) was then calculated and finally the total 
resistance {Rtot) = 1/(G) plotted against the computed 
conductance Gtun of the isolated tunnel junction. Since 
the average conductance of the diffusive region {Gdiff) 
is also known, equation (|l|) can be evaluated to yield the 
corresponding analytical result, curve (a) of figure ^. 

This demonstrates that in the range of validity of 
quasi-classical theory, quantitative agreement with the 
numerical scattering approach is obtained. 

V. NUMERICAL RESULTS FOR ANDREEV 
INTERFEROMETERS . 

Having examined a simple N-I-S structure, we now 
compare numerical results for the tight-binding structure 
of figure ||b, with the predictions of quasi-classical theory 
for the one-dimensional system of figure |^(a). The latter 
comprises a tunnel junction connected by diffusive 1-D 
wires to a fork. Each of the two arms of the fork is a dif- 
fusive wire, connected via tunnel junctions to infinitely 
long superconductors. The conductance of the diffusive 
wires is assumed to be much greater than that of the 
tunnel junctions. 

The two-dimensional tight-binding realisation of this 
structure is shown in figure ^(b) , which consists of a tun- 
nel junction (1) lying next to a diffusive region which is 
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in turn adjacent to two superconductors. The supercon- 
ductors are separated from each other by an insulating 
layer and from the diffusive region by two identical tun- 
nel junctions (2). The superconductors i = 1,2 have or- 
der parameter phases (jji , but are identical in every other 
respect. In order that they may successfully represent 
superconducting probes of infinite length, they are cho- 
sen in such a way that quasiparticle transmission through 
them is negligible. 

The diffusive region is of length Ldifj sites, each tunnel 
junction is Ltun sites long and the superconductors have 
a length Lsup- To model a superconducting reservoir, 
Lsup is again chosen sufficiently large such that there is 
negligible transmission through the superconductor. The 
system width and the width of both the diffusive region 
and the tunnel junction (1) is M sites. On the right of 
the diffusive region, the three insulating layers are each 
one site thick and therefore the superconductors are each 
of width M' (where 2M' = M - 3). 

In the simulation, the conductances Gi and G2 of the 
tunnel junctions are fixed at values which replicate the 
three situations of figure |4[ namely Gi = G2, Gi ^ G2 
and Gi ^ G2 to enable comparisons to be made with 
the analytic results. In each case, the phase difference 
between the two superconductors is varied and the to- 
tal conductance G plotted as a function of phase for 
a particular realisation of disorder in the diffusive re- 
gion. Ensemble-averaging over many disorder realisa- 
tions yields the conductance (G) which is independent 
of the microscopic configuration of the system and may 
be usefully compared with the results of equation (^) (see 
below) and figure ^. 

In what follows, we examine a sample of width M = 15, 
with a diagonal matrix element eo = 0.2 and periodic 
boundary conditions, for which the number of open chan- 
nels is iV = 13. The disorder is chosen to be = 1 and 
again from a graph of the form of figure 5, we obtain 
a mean free path of I « 4.9. In most cases of experi- 
mental interest, the conductance of the diffusive 'wires' 
may be considered to be much greater than that of the 
tunnel junctions. In order to take into account this situ- 
ation in our numerical simulations, one has to put some 
restrictions on the length of the diffusive region, since 
the conductance decreases with length. A compromise 
must therefore be found between the desire to increase 
the length into the diffusive regime and the wish to de- 
crease it in order to maintain a high conductance. In 
what follows a choice L^iff = 10 is made, for which 
(Gd^ff) « 5.1. 

To create a tunnel barrier of length Ltun = 1 and 
width AI, all diagonal elements Ci of sites within the 
barrier were set to a value = eg + Cfc. For an iso- 
lated tunnel junction (1) of width M — 15, the val- 
ues Cfc = 3.52, 12.27 yield respectively the conductances 
Gi = 2.0, 0.2 and for an isolated tunnel junction (2) of 
width M' = 6, the values = 1.87, 7.75 yield conduc- 



tances G2 = 2.0, 0.2. These values of Cb were used in 
the simulations of figure 9. As in the previous section, 
the choice eo = 0.2, Aq = 0.05, Lsup — 100 was made. 
For the superconductor alone, connected to crystalline 
normal leads of width M' = 6, the values of the normal 
and Andreev reflection and transmission coefficients were 
found to be: Rq = 0.02155, Ra = 4.92905, To = 0.04937, 
Ta — 0.00003. For such a structure, there are 5 open 
channels and as a consequence, the sum of these four co- 
efficients is 5. Finally, in order to model an insulating 
barrier, the diagonal matrix elements referring to a 
site i on the barrier between the superconductors, were 
each set to the large number £i = eo + 50. 

By combining the above components to yield the com- 
plete structure of figure ||(b) , one obtains the structure to 
be analysed numerically, whose parameters are as follows: 
total width M = 15, superconductor width M' — 6, num- 
ber of open channels in normal lead iV = 13, band filling 
eo = 0.2, chemical potential ^ = 3.8, length of tunnel 
junctions Ltun = 1, length of diffusive region Ldiff = 10, 
disorder width W ^ 1.0, conductance of diffusive region 
(Gdiff) = 5.1, length of superconductor Lsup = 100, su- 
perconducting coherence length ^ = 76 , superconducting 
order parameter Ao = 0.05. 

To carry out the simulation, the conductances Gi of 
tunnel junction (1) and G2 of tunnel junction (2) were 
fixed and a particular realisation {ei} of disorder in the 
diffusive region was selected. Then the phase difference 
(j) ~ — (t>2 between superconductors 1 and 2 was varied 
from zero to 27r. This was done by fixing 01 = and 
choosing 50 evenly spaced values of (j>2- For each value 
of 0, 200 different diffusive regions were obtained and 
the conductance G of the whole system computed for 
each. The ensemble-averaged conductance (G) was then 
calculated. The graphs of figure || show plots of (G) as a 
function of the phase difference (j> for the following four 
combinations of Gi and G2. They are (a) Gi = 0.2 
G2 = 0.2, (b) Gi = 2.0 G2 = 0.2, (c) Gi = 0.2 G2 = 2.0, 
(d) Gi = 2.0 G2 = 2.0. 

Apart from the different vertical scales, the numerical 
results of figure ^ and the analytic results of figure ^ share 
many qualitative features and also exhibit some interest- 
ing differences. Figure H(d) is comparable with ^(d); each 
exhibits a zero-phase minimum and a further minimum at 
(f> = TT. Similarly |9|(b) is comparable with ^(b); each ex- 
hibits a zero-phase maximum, with a minimum at = tt. 
The remaining curves compare less favourably. Whereas 
the analytic results of figures ^a) and ^d) are necessar- 
ily identical, there is no such restriction on the numerics 
and as shown in figure |^(a), decreasing the conductances 
Gi and G2 can produce qualitative changes. As a con- 
sequence, figure |^(a) possesses a zero-phase maximum, 
whereas figure Q(a) possesses a zero-phase minimum. 

Figures ^(c) and ^(c) also reveal some differences. 
Each possesses a zero-phase minimum, but at = tt, 
where the analytic result vanishes, the numerical result is 
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almost maximal. The inset in figure ||(c) shows a 'blowing 
up' of the region tt - 0.15 < < tt + 0.15, with d = 0.2 
and G2 = 2.0. The inset shows three curves, obtained by 
averaging over different numbers of disordered samples, 
namely 200, 1000 and 2000 realisations of the disorder. 
These demonstrate that in contrast with eq.(||), the nu- 
merical results ^(c) possess a shallow, local minimum at 

(j) = TT. 

Finally, we end this discussion by noting that for sys- 
tems with a small number of open channels, the be- 
haviour of an individual sample can be very different from 
that of the mean. For each of the four cases (a) to (d), fig- 
ure |o| shows each of the 200 plots of conductance G from 
which the ensemble averages of figure ^ were calculated. 
Apart from the case 10(b), where individual members of 
the ensemble behave in the same manner as the ensemble 
average, the nature of the extrema at (/) = 0, tt depends 
on the miscroscopic realisation of the disorder. We also 
note that by changing the dimensions of the sample, one 
can change the details of figure 9, but not the qualita- 
tive shape. For example by increasing the length Ldi// 
from 10 to 15, the local minimum at (/) = tt in figure 9c 
becomes more pronounced, but further increasing Ldiff 
to 20 causes the minimum to become more shallow. 

Having compared the quasi-classical theory of refer- 
ences 0- JT^ with the numerical scattering approach in 
the diffusive limit, we now examine the cross-over to the 
ballistic regime, where the former is inapplicable. We 
focus attention on the interferometer of figure 3(b) and 
examine the change in behaviour as the length Ldiff of 
the diffusive region becomes smaller than the elastic scat- 
tering length I. Apart from the change in Ldiff all other 
parameters are fixed to the values used in figure 9. Figure 
1 1 shows results for a diffusive region of length Ldi ff = 5 
and figure 12 for a length Ldiff = 1- Remarkably, apart 
from the overall increase in the conductance, the quali- 
tative shape of the curves is unchanged, despite the fact 
that the restrictions on quasi-classical theory are violated 
by these structures. 



is a quasi-one-dimensional theory and for the symmet- 
ric structures of figure 3, necessarily predicts a vanishing 
conductance at = tt. This symmetry is not present 
at a microscopic level and therefore there is no such re- 
striction on the results from an exact solution of the Bo- 
goliubov - de Gennes equation. Figures ^(a,b,d) suggest 
that for certain structures, this microscopic symmetry- 
breaking may be unimportant, but for other strengths 
of the tunnel barriers, figure ^(c) suggests that this ar- 
tifact of quasi-classical theory will not be observed ex- 
perimentally. In the clean limit, figures 11 and 12 show 
that although the overall conductance is increased, the 
qualitative shape of the conductance-phase curves is un- 
changed. 

The above results demonstrate that quasi-classical the- 
ory yields the correct shape for the ensemble-averaged 
conductance even down to extremely small system sizes 
and that results obtained for dirty systems can be appli- 
cable, in some cases, even in the clean limit. They also 
demonstrate that even without attempting a systematic 
extrapolation to a large number of channels, numerical 
multiple-scattering calculations on small structures, can 
yield results for ensemble-averaged properties of much 
larger systems. 

The results of figure |lj demonstrate that the essential 
properties of these interferometers are unchanged in the 
clean limit and therefore as already noted in [Q disorder 
is not a necessary feature of large amplitude Andreev 
interferometers. 
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VI. DISCUSSION 



In this paper, for the first time, a detailed comparison 
between quasi-classical theory and numerical multiple- 
scattering calculations has been carried out. To ensure 
that the simulated structures fall within the parameter 
range where the approximations of quasi-classical Green 
function methods hold, we have painstakingly examined 
each component of a given structure. For the N-I-S struc- 
tures of figure 1, figure 2 shows that there is quantitative 
agreement between the two methods. For the interfer- 
ometers of figure 3, there is broad qualitative agreement, 
although as shown in figures ^ and ^, some interesting 
differences are present. The theory of references m- |12] 
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FIG. 1. (a) Schematic picture of an N-I-S junction for anal- 
ysis using circuit theory, (b) Picture of the tight-binding lat- 
tice used to model a two-dimensional N-I-S structure. 

FIG. 2. Conductance of the N-I-S structures of figure 1, 
as a function of the conductance Gtun of the tunnel junc- 
tion. The various curves refer to: a) analytic theory; b) 
numerical simulation; c) resistance of the diffusive region 
< Rdiff >= 0.63 =< Gdiff >^^; d) asymptotics at high 
transparency Gtun ^< Gdiff >', e) asymptotics at low trans- 
parency Gtun <:< Gdiff >■ 

FIG. 3. (a) Schematic picture of an interferometer an- 
alyzed using circuit theory, (b) Picture an interferometer 
formed from a two-dimensional tight-binding lattice. 

FIG. 4. Analytic results for the conductance versus phase 
in the diffusive Umit: a) Gi = 0.2 and G2 = 0.2. b) Gi = 2.0 
and G2 = 0.2. c) Gi = 0.2 and G2 = 2.0. d) Gi = 2.0 and 
G2 = 2.0. 

FIG. 5. < Gdiff > Laiff/N versus L^iff. The plateau 
region signifies the diffusive regime. Here the number of open 
channels is A'' = 9. In the inset standard deviation is also 
shown. 

FIG. 6. Conductance Gtun of the tunnel junction as a 
function of the barrier height et- 

FIG. 7. The conductance G(Ao) of a clean superconduct- 
ing region of length Lsup = 5, width A4 — 50, plotted as a 
function of the site energy eo, for 5 different values of Aq. 

FIG. 8. Plots of the mean conductance 

(G(A())) — N — Ro+Ra as a function of Aq, for a diffusive con- 
ductor of width M = 10 and four difi^erent lengths. The left 
insert shows the conductance < g >= (G(A()))/(G(0)), scaled 
by the normal state conductance. The right insert shows the 
rms deviation a = ((G(Ao) - (G(Ao)))^)^/^ 

FIG. 9. Numerical results for the conductance versus phase 
for Ldiff = 10: a) Gi = 0.2 and G2 = 0.2. b) Gi = 2.0 and 
G2 = 0.2. c) Gi = 0.2 and G2 = 2.0. d) Gi = 2.0 and 
G2 = 2.0 
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FIG. 10. Numerical results for the conductance versus 
phase for individual realizations of the disorder : a) Gi = 0.2 
and G2 = 0.2. b) Gi = 2.0 and G2 = 0.2. c) Gi = 0.2 and 
G2 = 2.0. d) Gi = 2.0 and G2 = 2.0. 



FIG. 11. Numerical results for the conductance versus 

phase for L^,// = 5: a) Gi = 0.2 and G2 = 0.2. b) Gi = 2.0 
and G2 = 0.2. c) Gi = 0.2 and G2 = 2.0. d) Gi = 2.0 and 
G2 = 2.0. 



FIG. 12. Numerical results for the conductance versus 
phase for Luff = 1: a) Gi = 0.2 and G2 = 0.2. b) Gi = 2.0 
and G2 = 0.2. c) Gi = 0.2 and G2 = 2.0. d) Gi = 2.0 and 
G2 = 2.0. 
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